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We calculate the scaling exponents of the two-dimensional correlated percolation cluster's hull 
and unscreened perimeter. Correlations are introduced through an underlying correlated random 
potential, which is used to define the state of bonds of a two-dimensional bond percolation model. 
Monte-Carlo simulations are run and the values of the scaling exponents are determined as functions 
of the Hurst exponent H in the range —0.75 < H < 1. The results confirm the conjectures of earlier 
studies. 



I. INTRODUCTION 

The world around us is chaotic and seemingly random, 
but one can still find regularities. Take mountain ranges 
— these irregular jagged structures may seem intractable 
to analysis, but often display an interesting property — 
they look the same at different length scales — they are 
self-similar. This phenomenon is not limited to the sur- 
face features of planets, but is also found in many other 
places — deposited metal films ripple- wave turbu- 
lence 13], crack fronts in material science 0-0] > cloud 
perimeters Q, passive tracers in two-dimensional fluid 
flows [1-GJ], etc. 

The analysis of these physical systems is often reduced 
to determining the scaling exponents characterizing the 
rough surfaces involved. For example, the convective- 
diffusive transport of a passive scalar in a random two- 
dimensional steady flow is determined by the scaling ex- 
ponent of the isolines of the underlying stream-function 

The aim of this paper is to numerically calculate the 
values of these exponents depending on the roughness pa- 
rameter H . In Sec. [Ill we start off by giving an overview 
of the concepts used — rough surfaces and percolation 




Rough surfaces (H > 0) 



Correlated surfaces (H < 0) 



Figure 1: The potential (|4|) is made up of components with 
different amplitudes. The amplitude of a component of scale 
A is proportional to X H . For H > the wider "hills" start 
to dominate the landscape and once H > 1 only the widest 
"hill" matters. Conversely, for H < local fluctuations start 
to gain in influence and once H < —0.75 the wider hills lose 
any influence on the scaling exponents of the surace. 
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Figure 2: An example of a random potential (H = and so 
Dh = f.5 [3,@|) with isolines separating the area into "land" 
and "sea". 



clusters, what we mean by correlations, and a mapping 
between the two classes of models. Numeric calculations 
are done through the Monte-Carlo simulations using the 
two-dimensional bond percolation model; the procedure 
is described in Sec. Mil Interpretation of the resulting 
data requires overcoming the finite size effects and the 
convergence problems, which is the subject of Sec. IIVI 
Finally, Sec. [V] provides a brief summary of the results 
and a future outlook for the studies of the statistical to- 
pography of random surfaces. 



II. OVERVIEW 

Rough surfaces. Let ip(x,y) = ip(r) be the height or 
potential of a self-similar random two-dimensional sur- 
face. We define the roughness exponent H — also known 
as the Hurst exponent — through the surface height drop 
at distance a = |a| fl3j : 



[V>(r) - tp(r + a)] ) cx |a| 



t 2H 



(1) 



where angular braces denote averaging over the coordi- 
nate r (or also over an ensemble of surfaces). This scaling 
law assumes that ao < |a| < oi, where ao and a\ are the 
lower and upper cut-off scales, and < H < 1. Relation 
([!]) also describes self-similarity — the height of a "hill" 
on the surface is a power law of its diameter, so hills at 
different scales have the same proportions. 



A more generic description, which is not limited to the 
positive values of H, can be given through the power 
spectrum Pk We assume that 

(Vk)=0, (^Mv)=Pk<W, (2) 
and define the spectral density Pk as a power law: 

P k oc|kr 2ff -\ for |k|<a^. (3) 

This also allows us to conveniently divide the "multiscale" 
potential ip(r) into a sum of "monoscale" functions, 

1 p(r)= M*), (4) 

A i =2 i A 

where each component in the sum represents a function 
with a single characteristic length. The effect of the pa- 
rameter H on the amplitude of the "monoscale" compo- 
nents of the "multiscale" potential can be seen in Fig. [T] 
The potential ip can be also characterized through a 
correlation function (covariance) 

C(a) = (V(x))(V(x + a)). (5) 

Indeed, for potentials conforming to (J3)) , it is a power law 

(7(a) cx \sl\ 2H a Q < |a| < oi, (6) 

which is valid for the range — 3/4< H < 1. 

With a specific height h, the random potential defines 
a set of isolines ip(r) = h (see Fig.[2]for an example). The 
height h can be interpreted as the "sea level". So, a single 
isoline can be looked at as the coastline of an island or 
a lake. The coastline is also a self-similar structure, and 
one can quickly run into difficulties when trying to mea- 
sure its length — the result depends on the size of the 
measuring stick used (l5| . The parameter that best char- 
acterizes isolines is their scaling exponent — also called 
their fractal dimension Dh — which is a fractional num- 
ber. The length of a coastline can then be given through 

(iWlocA 1 ^, (7) 

where A is the size of the measuring stick. Dh is a nontriv- 
ial function of the underlying surface's Hurst exponent: 

D h = D h (H). (8) 

Alternatively, we can take two points at a distance a on 
a coastline, and calculate the length of the line between 
them, assuming a fixed measuring stick A = Ao: 

(L(a))(xa Dh . (9) 

The percolation problem is concerned with the 
structures that form by randomly placing elementary ge- 
ometrical objects (spheres, sticks, sites, bonds, etc.) ei- 
ther freely into continuum, or into a fixed lattice (Fig. [3]). 
Two objects are said to communicate, if their distance 
is less than some given Ao, and communicating objects 
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Figure 3: Examples of regular percolation lattices. Square 
lattice site model to the left and square lattice bond model to 
the right. Largest clusters have been circled. 




Figure 4: An example square bond percolation cluster. The 
cluster is made up of the bold segments, the zig-zag is the 
hull and the dashed line is the unscreened perimeter. 



form bigger structures called clusters. Percolation the- 
ory studies the formation of clusters and their proper- 
ties. The more interesting aspect is when and how is an 
infinite cluster formed. This depends on the lattice site 
occupation probability rj. The minimum site occupation 
probability when an infinite cluster appears is called the 
percolation threshold r/ c . Near this probability, the per- 
colation model displays a critical behavior and long-range 
correlations. 

Percolation theory is used to study and model a wide 
variety of phenomena: from a fluid flow in a porous 
medium to thermal phase transitions and critical behav- 
ior in magnetism with dilute Ising models. 

The structures that can be identified in conjunction 
with a percolation cluster are the cluster itself, the hull 
and the unscreened perimeter (Fig. [4} . Aside from these 
one more structure can be identified — the oceanic coast- 
line [l6||. It is comprised of the hull of a cluster and 
also the hulls of all the holes (reverse clusters) within 
it. Near the percolation threshold, all of these structures 
(the cluster, the hull, the unscreened perimeter and the 
oceanic coastline) are fractals and can be characterized 
by scaling exponents. 

Looking at a percolation cluster and an isoline of a ran- 
dom potential, we can identify similar structures: one can 
look at the percolation cluster as an island and its hull as 
the coastline of the island. This idea will be fleshed out 
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in more precise terms farther below, where a mapping 
between the two models is described. 

One very important and useful aspect about the scaling 
exponents is a phenomenon known as universality [14 1 - 
within specific universality classes, the scaling exponents 
take the same value across different percolation models. 
More specifically, they are invariant to small fluctuations 
or distortions in the lattice structure (for instance, de- 
caying exponentially to distance). This means that the 
scaling exponents for both the random square bond lat- 
tice, and the random square site lattice are the same: the 
both models belong to the same universality class of the 
two-dimensional uncorrelated percolation. 

Correlated percolation. Percolation lattice does 
not need to be completely random, but can entail cer- 
tain correlations. Here we describe the percolation lat- 
tice through an infinite set of random variables 8i which 
are unity at occupied sites and zero at empty sites (i de- 
notes the site number). Then we can characterize the 
correlations through the correlation function 



P 



(10) 



where p = (9i) is the site occupation probability. 

Alternatively [I?} , correlations can be brought into the 
percolation model by assigning each lattice site a random 
number pi £ [0, 1] where (pi) = p. The site values are 
then calculated as 



Q(pi - Xi), 



(11) 



where O(x) is the Heaviside step function and {xt} are 
independent random variables uniformly distributed in 
[0, ll. The correlation function is 



c p (x t -*j)= (jPiPj) -p 2 



(12) 



and cg(a) = c p (a) = c(a) fl4l |. 

We are interested in algebraically decaying correlations 
so that 



c(a) ex 



,2H 



H < 0. 



(13) 



It is believed [l4| that the universality class is determined 
by the two-point correlation function. One can show that 
at H < —3/4, the model belongs to the universality class 
of uncorrelated percolation [17| . However, in the range 
—3/4 < H < 0, the correlations do affect the scaling 
exponents. 

There exists a simple mapping between the 
rough surfaces and the percolation model fill. Il9j. 

According to it, the local maxima of the potential define 
the lattice sites and the lattice bonds are obtained by 
drawing fastest-ascent paths from all the saddle points 
(Fig. [5]). The surface ^(x) is "flooded" at a given level h 
and a bond i is left connecting if the saddle point Xj on it 
is above the water (is land), that is when > h. As a 

result, we get an irregular two-dimensional lattice; recall 
that as per universality small distortions of the lattice 
don't affect the resulting scaling exponents. With this 




Figure 5: Mapping between the random surfaces and the per- 
colation problem. The bond between two maxima mi and 
?7i2 is connecting, if the saddle point S12 on them is above 
the flood level. So, saddle points are mapped into percolation 
bonds. 
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Figure 6: Known results, conjectures and interpolations for 
the scaling exponents as functions of the Hurst exponent H. 
While these functions can be approximated as linear func- 
tions, they are in fact non-trivial. 



mapping, we can relate the islands formed at flooding to 
the resulting clusters, and their coastlines to the hulls of 
the said clusters. Also, if the surface correlation function 
§5§ is a power law © , then so is the correlation function 
(JTUJ) for the percolation model (p~5|) , where the parameter 
H is the same. Due to universality, the scaling exponents 
of the matching structures are also the same. 

The scaling exponents are of interest in many applica- 
tions so, there is a need to calculate their values depend- 
ing on the underlying surface's roughness parameter H. 
In Fig. [5J we can see the known results (numeric and 
analytic) , and also interpolations and conjectures for the 
range —3/4 < H < 0. Our next task is to run simula- 
tions to numerically shed light on these gray areas (the 
scaling exponents of the hull and the unscreened perime- 
ter). As these exponents behave the same way for the 
both problems of rough surfaces and correlated percola- 
tion, we can calculate them using the model which is the 
most convenient from the numerical point of view. 



III. MONTE-CARLO SIMULATIONS 

Generation of percolation clusters. We calculate 
the scaling exponents using the model of correlated two- 
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Figure 7: For the i-th layer (r) in 0, we divide space into 
a grid of cells with side length of Ai = 2 1 (Ao = 1)- At the 
center of each cell j, a cone is placed with height ctj = rj2 lH , 
where rj is an independent uniform random variable in the 
range [—0.5,0.5]. The cones have a diameter of 3A^ and so 
overlap. 

dimensional bond percolation on square lattices. Our 
first task is to generate said percolation models so that 
they conform to the correlation function (fTB")) . For this, 
we first generate random potentials with the requested 
roughness H . We take the "flood level" as the value of the 
potential at some starting location Xo, so that h — ^(xq), 
and use it to map the random surface model into one of 
the percolation models. This is done by overlaying the 
percolation lattice on the rough surface and calculating 
the bond values using the Heaviside step function as 

Ot = 6 (V>(xi) - V(x )) . (14) 

This approach is similar to the one described by equation 
(fTT|) and preserves the correlation exponent. It is slightly 
different from the maxima-saddle point mapping, but due 
to universality, there will be no change in the values of the 
scaling exponents. Using the potential as the underlying 
model allows us to get the results in the parameter range 
-0.75 < H < 1. 

To generate random potentials, we exploit formula 
((4]) . We generate different components (layers) of the po- 
tential for different lengths and sum them up. The «-th 
layer ^Ai( r ) is formed by a grid of cells, where the grid 
cell side length is 2 l , and at the center of each cell j we 
place a cone of height rj2 lH , where rj is an independent 
uniform random variable in the range [—0.5, 0.5] (Fig. [7]). 
The cones have a diameter of 3Ai, so the resulting over- 
lap yields a smoother potential. However, the overlap or 
the shape used (cone, in this case) does not significantly 
affect how fast the results converge to the asymptotic 
power laws. So, one can choose a shape more optimized 
for the speed of numerical calculations (for instance, one 
could use a simple block or a cylinder). The number of 
layers required depends on the value of H. For example, 
at H — 0, we can work with i = . . . 30. 

Gathering data is a matter of generating percola- 
tion clusters of different sizes and tracing the structures 
of interest within these. For the hull, this is done in the 
following way. We constrain ourselves to a a x a box, 
and start to dynamically trace the hull from the center 
point x until it reaches any sides of the box (Fig. [5J. 
Hulls that make a full circle are discarded. To get ad- 
ditional data, we also trace backwards from the center, 



CO 




a 



Figure 8: Calculating the length of a hull for size a. We start 
from the center of a a x a box and dynamically calculate bond 
values as we trace the hull until reaching the side of the box. 
We discard hulls that make a full circle. The box sizes picked 
are a = 8, 16, 32, . . . , 1024, . . .. 




Figure 9: Counting the size of the hull can be done in many 
ways. One can simply count the number of segments (hull), 
the number of bonds the hull touches from both sides (sides), 
the total number of bounds the hull touches (bonds), the num- 
ber of times a single bond sticks out (ends), or the number of 
times the hull forms a straight line of 4 segments (lines). All 
of these scale with the same power law. 



so that the both ends of the hull reach the box's sides. 
We do this millions of times for differently sized boxes 
a = 8, 16, ... , 1024, . . .. While for some cases, the range 
of values a = 8, 16, . . . , 512 is sufficient, slow convergence 
in some regions of H forces us to use larger lattices. How- 
ever, the size is limited by computational resources and 
in our case it was not practical to go over a — 2048. 
Once we have a hull, we want to measure its size 
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(length). While the easiest way would be to just sum 
the number of segments in the trace line, it is also possi- 
ble to determine the size by other properties (Fig. |9|) . All 
of these properties scale with the same power law. We 
can denote the size of a single hull as Li(a,j), where i indi- 
cates the property. As this is different for each individual 
hull, we find the average value 



d{a ) = {Li(aj)) , 



(15) 



and as per equation (j9]) , this should scale as a power law 
with the exponent Dh- However, the scaling is asymp- 
totic (a — > oo), so for finite a, there are sizable deviations 
called finite size effects. These come from the geometry 
and finite size of the lattice. Using (j9j), we can estimate 
the exponent as 



ln 2 



C(a,j-i] 



(a.j = 2 a j. 



i)- (16) 



Plotting this for the uncorrelated percolation (Fig. [TIJ)l . 
we can see how the different properties converge towards 
the value Dh — 7/4. The finite size effects are strongly 
manifested for the smaller lattices. 



IV. DATA ANALYSIS 

To calculate the scaling exponents from the data, 
the following assumptions are made: 

1. the mathematical expectation for each property can 
be described as an infinite series 



Li(a) = ^A ltl a a ^', 
^=1 



< a 



(17) 



where i = 1, 



. , m; 



2. a Ul = a M for fx = 1, . . . , m; 

3. the leading terms in the sum are linearly indepen- 
dent (det IIA^II ^ 0). 

After these assumptions, we apply a variation of the least 
squares method described in [20( and previously used in 
[16j . The method works if the assumptions made are 
correct (the method also validates them) and yields us 
the value a\, which is the scaling exponent we are look- 
ing for. So, the reason why we counted all the different 
properties for the hull (Fig. [9]) , is that they are necessary 
for this method. 

For the unscreened perimeter, a similar approach is 
taken. The unscreened perimeter is obtained by taking 
a hull but pruning it from "fjords". 

Convergence problems. In some areas the calcula- 
tions are hindered by very slow convergence. Let param- 
eter S represent the length at which the covariance of the 
generated potential starts to differ from that of the ideal 
Ial : 



l2H law: 



bonds - 



Figure 10: Convergence of the hull properties for the uncor- 
related percolation (towards Dh = 7/4 = 1.75). 




S = 3 ■ 2 s 



(18) 



Figure 11: Covariance of a potential versus the power law 
x 2H . Here S marks the length at which the covariance starts 
to diverge from the power law. 



where s indicates the smallest-scale layer index (Fig. ITT]) . 
When S decreases (by adding bottom layers), local fluc- 
tuations start to gain in influence compared to those of 
the long-range correlations. This causes a strong finite- 
size effect and the scaling exponents behave as if H was 
smaller (Fig. [T2"|) . Conversely, when increasing S, the 
scaling exponents can initially behave as if H was greater 
than it really is. 

To get over these distortions, one could calculate for 
bigger lattice sizes. But often this is not an option as con- 
vergence can be very slow and computational resources 
are limited. Another way would be to manually find 
the optimal layer configuration that minimizes distor- 
tions. This is the approach we took and yielded good 
results for the hull (aside from the values H = —0.75 and 
H = +1.00). However, the convergence of the unscreened 
perimeter is very sensitive to changes in the layer con- 
figuration and for most data points did not yield clear 
results. 

The results can be seen in Fig. Q2] and Tab. [TJ The 
hull behaves as expected. While it did not yield clear 
results at H = —0.75 and H = 1.00, the extrapolations 
provided seem to indicate that it terminates at 7/4 and 
1 respectively. Due to the convergence problems, the re- 



6 



S-2.000 
S-4.500 I 
S=10.125 
S-22.781 



Figure 12: Convergence of a property of the unscreened 
perimeter at H = (towards theoretically known D u = 1.5) 
for different values of parameter S. 



suits for the unscreened perimeter are not as clear. How- 
ever, one can say that at least that the result do not 
contradict the analytical findings and support the ap- 
plicability of a nearly-linear interpolation between the 
points D u (0) = 1.5 and D u (-3/4) = 4/3. 



V. CONCLUSION 

We have run Monte-Carlo simulations to determine 
the scaling exponents of the hull and the unscreened 
perimeter as functions of the Hurst exponent in the range 
—0.75 < H < 1. For this, we first generated random po- 
tentials conforming to the required correlation function 
by summation of component potentials of different char- 
acteristic lengths and mapping the potential into per- 
colation models. Hulls and unscreened perimeters were 
traced from these models, and their lengths calculated for 
different scales. A variation of the least squares method 
was used to obtain the values of the exponents. 

The results confirm the previously known data in the 



range < H < 1 and also the conjectures for the behav- 
ior in the range —0.75 < H < 0, see Fig. [13] The par- 
ticular results regarding the fractal dimension Dh(H) of 
hulls confirm that for < H < 1, the 4- vertex model (i.e. 
the rough surfaces in 1+1-dimensional geometry) [2l| be- 
longs to the same universality class as the isotropic Gaus- 
sian self-affine surfaces (assuming the respective equality 
of the Hurst exponents). Indeed, comparing the numer- 
ical results of Ref. |2l| and those of the current study 
shows that for the entire range of < H < 1, the values 
of Dh(H) coincide within the uncertainties of ca 10~ 3 . 
An important consequence is that the conjecture about 
the super-universality of the loop correlation exponent 
xi(H) = i 01 3 ( which has been exploited in several 
studies, c.f. |22H25j |) is clearly rejected: this conjecture 
implies Dh{H) = | — which falls well beyond the un- 
certainty margins of the present simulation results (for 
instance, at H = ^, the conjectured value | falls far 
from the range of 1.2820 ± 0.0008). 

The obtained results are valuable for a range of practi- 
cal applications (such as the turbulent transport in quasi- 
stationary velocity fields), for which the scaling expo- 
nents have been analytically expressed via the fractal di- 
mension of the hull. As a future outlook, our method can 
be applied to calculate other scaling exponents of the cor- 
related percolation problem and statistical topography, 
such as the fractal dimensions of the clusters (oceanic 
coastlines), percolation backbone, etc. It can be also 
extended to study the scaling laws of the transport on 
quasi-stationary velocity fields. 
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